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I. INTRODUCTION 

The recent discovery [HI^EiS,01 °f the + baryon with a mass of 1540 MeV close to the values predicted in Refs. 
0, has led to considerable excitement in hadronic physics. One of the main ingredients involved in the prediction 
of G + 00 is the limit of the large number of the QCD colors N c and the l/iV c -expansion. In spite of the numerical 
agreement between the theoreti cal p rediction of Refs. |(J, an d the experimental data, recently certain doubts have 
been expressed in Refs. [8ll^. ll0llll| concerning the consistency of the methods used in Refs. 0,0 with the systematic 
1 /iVc-expansion. 

As is well known, the standard picture of large- iV c baryons assumes that they can be described by a mean-held 
solution of some effective theory equivalent to the large- N c QCD 12] . The precise action of this effective theory 
is not known. Nevertheless, the general properties of large- N c baryons can be derived from the symmetry of the 
mean-field solution. It is assumed that the mean-field solution violates both the flavor SU (3) group and the group 
SO(3) ~ SU(2) of usual space rotations. However, the mean field is supposed to be invariant under combined SU(2) 
rotations in the flavor and in the usual space. The realization of this symmetry can be found in various models ranging 
from the naive quark model to the Skyrme model 0, 0, |^. At large N c these models lead to certain relations, 
which are not sensitive to peculiarities of different models. These model-independent relations are believed to hold 
also in the large- N c QCD. Many of these results can be derived using large- N c consistency relations [ltl ITtI ITsl ITsj. 

In the approaches based on the mean-field solutions of effective theories, one should properly treat the symmetries 
violated by the mean field. In principle this problem can be solved in every order of the l/iV c -expansion. However, 
sometimes the systematic large- N c analysis is replaced by the rigid-rotator approximation. In the chiral limit the 
rigid-rotator Hamiltonian has the simple form 

Hmt = y A + y A. (i) 

a— 1 a— 4 

where J a are generators of the SU(3) group and If. are the "moments of inertia" of the rotator. Although the structure 
of the effective rigid-rotator Hamiltonian JQ) is inspired by the large- N c limit, the rigid-rotator approximation is not 
equivalent to the systematic l/A^-expansion. Therefore the results based on the rigid-rotator approximation should 
be taken with certain caution. 

The analysis of exotic states in Refs. 0,0 was based on the effective Hamiltonian JJJ. The moments of inertia Ik 
grow with N c : 

h = 0(N C ) . (2) 

Therefore the low- lying states described by the Hamiltonian (JTJ have the excitation energy 0(N~ 1 ). These low- lying 
0(N^ 1 ) excitations have nonexotic quantum numbers coinciding with the quantum numbers which can be obtained 
in the naive quark model with N c colors. 

In Ref. 20| it is argued that the rigid-rotator Hamiltonian ^ can be applied to describe exotic states whose 
excitation energy is 0(N®). On the other hand, the authors of Refs. pjl IsL fiol HH] claim that the extension of the 
rigid-rotator approach to the exotic sector is inconsistent with the systematic l/iV c -expansion. 
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The argumentation of Ref. |20| includes a toy model of exotic states (electrically charged quantum particle moving 
in the field of a static monopole) which was studied earlier in Ref. [21]. Within this toy model the relevance of the 
rigid-rotator approximation for the description of the exotic states was explicitly checked in Ref. 20] . On the other 
hand, the objections of Refs. [8LI9L IT0IITH are mainly based on the fact that the 0(A^°)-excitations should be described 
by the quadratic form of fluctuations around the mean field and not by the effective rotational Hamiltonian. 

In order to clarify this controversy, we have chosen several toy models which reproduce many properties of QCD, 
including the large- N c counting for baryons. In these models we compare the rigid-rotator approximation with the 
systematic 1 /A^-expansion. 

One of our models is exactly solvable. In contrast to the toy model considered in Ref. |2£j, the exactly solvable 
model studied in this paper reveals several inconsistencies of the rigid-rotator approximation. In short, these problems 
can be described as follows. 

1) The rigid-rotator approximation produces correct quantum numbers for the low-lying (nonexotic) 0(N~ 1 ) exci- 
tations but leads to wrong values of the energy. 

2) The spectrum of the rigid-rotator Hamiltonian contains 0(N®) exotic excitations which are absent in the exact 
solution of the toy model. 

One can also construct a systematic 1 /A^-expansion without using the exact solution. In this paper we compute the 
two first orders of the 1 /A^-expansion for the energy spectrum. We also find the spin-dependent part of the next-to- 
next-to-leading (NNL) order. Since our model contains fermion variables, the equations appearing in the systematic 
1 /iV c -expansion are quite similar to the well-known equations of the many-body physics [22| . In particular, the leading 
0(N C ) order of the l/A^-expansion is described by the Hartree equation. The next-to-leading (NL) 0(N®) order leads 
to the large- N c version of the RPA (random phase approximation) equations. The spin dependent part of the NNL 
0(N~ 1 ) order is described by the Thouless-Valatin formula [2^, whereas the rigid-rotator approximation leads to 
the Inglis formula 24] for the same quantity, which shows the incompatibility of the rigid-rotator approximation with 
the 1/TVc-expansion. 

The main part of the paper is devoted to the models containing colored fermions. In addition we consider another 
class of models based on local actions with color-singlet boson fields. We show that in a general case the rigid-rotator 
approximation does not work in the exotic sector of these models either. 

The structure of the paper is as follows. In Section [HI a brief description of the spectrum of the rigid-rotator 
Hamiltonian is given. In Section lllll we describe the general structure of the l/A^-expansion for models with colored 
fermions in terms of Hartree and RPA equations. In Section IIVI we introduce an exactly solvable toy model with 
colored fermions. The exact solution of this model is described in Section In Section fVTI we construct the 1/N C - 
expansion for the energy spectrum of the toy model in the leading and next-to-leading order. The spectrum in the 
NNL order is studied in Section IVIII The discrepancy between the rigid-rotator model and the systematic 1 /N c - 
expansion is discussed in Section lYlIII In Section ITxl we revisit the toy model of Refs. H3,|^] where the rigid-rotator 
approximation was successful, we show that this success is accidental. In Conclusions we summarize the results 
obtained in the toy models and make comments on the relevance of these results for the exotic baryons in large- N c 
QCD. 



II. RIGID-ROTATOR APPROXIMATION 



We begin with a brief review of the properties of the rigid-rotator Hamiltonian QJ. This Hamiltonian acts in 
the space of the wave functions ip(R) depending on the SU(3) matrices R. Operators J a appearing in Eq. are 
generators of the right SU (3) rotations 

[exp (iu a J a ) ip] (R) =^[R exp(iw A a /2)] (a = 1,2,..., 8). (3) 

Here A a are Cell-Mann matrices. The components Ji.,2,3 are directly interpreted as spin operators, whereas the 
spectrum of the component J& is constrained by the condition 14] 

j8 = -^r (4) 

where B is the baryon charge. 

The left rotations with generators T a 

[exp {iu a T a ) xj}} (R) = if, [cxp(-iuj a X a /2)R} (5) 

are interpreted as SU(3)q flavor transformations. 
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In the space of the rigid-rotator wave functions one nas the identity 



Vtt - Vj j 



(6) 



In terms of the standard (p, g) nomenclature of the irreducible SU(3) representations, generators T a and J a belong 
to complex conjugate representations (p,q) and (q,p) respectively. The Casimir operator is 



£ T a T a = J2 Ja Ja = C S 2 U ^ (p, q) = C S 2 U ^ (q, p) = \ (p 2 + g 2 + pq) + p + q. 



(7) 



Together with the Jg constraint (@J this leads to the well-known problem of the 1 /^-expansion for the baryonic 
states in the case of the SU(3)a flavor group: in the large- iV,- limit all traditional (octet, decuplet,...) and exotic ( e.g . 
antidecuplet) representations of SU(3)^ are lost. Instead, one has to work with the large- N c (p, q) representations |2q 
imitating octet, decuplet, antidecuplet etc. The condition 0} leads to the following constraint on the allowed (p, q) 
representations 



P + 2q > N c 



(8) 



In particular, the large- N c analogs of the main SU (3)a multiplets are 

1) octet: p — 1, q= (N c — l)/2, 

2) decuplet: p = 3, q = (N c - 3)/2, 

3) antidecuplet: p = 0, q = (N c + 3)/2. 

The eigenvalues of the effective Hamiltonian H Iot corresponding to the SU (3)a representation (p, q) and to the spin 
J (for B = 1 states) are 



£ ro t = 777- [p 2 + q 2 +pq + 3(p + g)] + 

DJ2 



1 1 

2h~2h 



J(J+l) 



N 2 

1 _ c 

24/ 2 



(9) 



Following Ref . |20j| , one can define the "exoticness" parameter (in order to avoid the confusion with the energy we 
replace notation E, used in Ref. |20j. by £) 



e= -(p + 2q-N c 



running over the nonnegative integer values £ = 0,1,2, 
Then 



E IC 



Let us consider the limit 



2/ 2 



■« + '> + 4(H + §(§+' 



i i 



J(J + 1). 



The spin J is constrained by the condition 



iV c — ► oo , p, £ = const . 



]£zfl<j<£±* 



In the case ( = flwe have 



so that the spectrum (|llfl has the form 



p + 2g = JV c , J=| (£ = 0) 



^ot = || + ^-J(J + l) (£ = 0). 
In the case of arbitrary £ we find from Eq. i|ll[l in the limit l|12l) 

s rot = ^(e + i) + o(iv- 1 ). 

This expression shows that the low- lying states have £ = 0. According to Eq. (|15|l the splitting between the £ = 
states is controlled by the J(J + 1) term and has the order 0(N~ 1 ). On the contrary, the states with £ ^ have 
0(N®) excitation energies. 



(10) 

(11) 
(12) 

(13) 
(14) 

(15) 
(16) 
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III. l/iV c -EXPANSION FOR FERMION SYSTEMS IN TERMS OF HARTREE AND RPA EQUATIONS 

In this section we describe the general structure of the 1 /iV c -expansion for fermionic systems with the Hamiltonian 



c A 1 A 2 A 3 A i \c' = l J \c=l / 



(17) 



Here -0,4c are fermionic operators ("quarks") with the "color" index c — 1,2, . . . , N c . Index A covers other degrees of 
freedom. We assume the standard fermion commutation relations 

{ipAc, V'l'c'l = ^AA'Sccf ■ (18) 

The Hilbert space of states can be constructed by acting with the operators i\)\ c on the "vacuum" |0) obeying the 
condition 

^a c |0)=0. (19) 

The hermiticity of the Hamiltonian implies that 

V A 1 A 2 \A 3 A i = VAiA 3 \A 2 A 1 - (20) 

We also assume 

^A 1 A 2 \A 3 A i = V A3 A 4 |A 1 A 2 • (21) 

The interaction 117|) is SU(N C ) invariant. We have no local gauge symmetry since the model has no space: the 
fermions are "pinned" at one point. The factor of 1/N C is separated in the RHS of Eq. (| 1 T|l in order to keep V fixed 
in the limit N c — > oo. 

One way to construct the l/_/V c -expansion is to use the traditional operator approach. Then one arrives at the 
large- N c version of the Hartree and RPA equations. These equations are well known in the many-body physics 
as well as their connection with the 1/A-expansion j2||. Alternatively one can derive the 1 /Ac-expansion using the 
path integral approach which is described in Appendix lAl 

At large N c the energy E of the low-lying states can be expanded as 

E = N c Ei +E + N~ 1 E- 1 + ... (22) 
In the leading order of the 1 /iV c -expansion the ground state energy is given by the Hartree approximation: 

Ei = \v ABWB ,P { ° c A c) P { °$. (23) 

Here P B °^ is the projector onto the single-particle occupied states 

P ba C) = £ 4> { B<t>A* ' (24) 

i: occ 

which can be found by solving Hartree equations: 

TTAB^ = £«4 n) . (25) 



kab = Vab\a>b> r*<>B- = Vab\a>b>P b °7 ■ ( 2 6) 



?: occ 



In order to find the 0(N®) correction to the energy of the ground state and the energy of the 0(N®) excitations 
one has to solve the large- N c version of RPA equations 

y ] y , (A m injX n: j + B m i n jY n j^j = Q u X mi , (27) 

n:nonocc j:occ 
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n:nonocc j:occ 

Here indices i,j correspond to occupied Hartree states whereas indices m, n stand for nonoccupied states. We use the 
notation 

-Bminj = ^Vminj j (^0) 

W a6od = ^ ^^bW^Vab^d- (31) 

Although equations (|27ll . (|28|) are quite similar to the standard RPA equations (see e.g. Ref. [I^]), Eqs. I|27l) . i|28|) 
have several different features. First, the color indices do not appear in Eqs. l(77l) . Second, we have dropped 

those terms of the standard RPA equations whose contribution is l/iV c -suppressed. Next, our RPA equations are 
formulated in terms of the Hartree and not Hartree-Fock occupation picture (the Fock term is a l/A^-correction 
which can be treated separately). The path integral derivation of Eqs. 1|27|) . 128|) can be found in Appendix lAl 

The RPA equations (|27fl . (|28[> can be written in the compact form 



where 



SZ" = n„Y,Z u , (32) 



^=(^), S = ( J _° x 1 , (33) 



Using Eqs. © and QQJ, we find 



Once the RPA equations H27|) . I|28l) are solved, the 0(N®) part i?o in the expansion ll'L'li of the energy equals 

E ° = -\ E E - e «) + E n - + ^ ■ ( 36 ) 



m:nonocc i:occ w ^ 

In the last term on the RHS the summation is over the positive RPA eigenfrequencies f2„ > (if any) . The harmonic 
excitations are parametrized by the integer numbers n v = 0,1,2,... The expression (|36H for Eq includes several 
contributions. In addition to the pure RPA term it contains the Fock contribution and the correction corresponding 
to the fact that our Hamiltonian l|17fl is not normally ordered. 

The calculation of the 0(N~ l ) contribution to the l/A^-expansion 122(1 requires more work. In this paper we 
consider only the spin-dependent part of the 0(N~ 1 ) contribution for so called rotational excitations. This part of 
the 0(N~ 1 ) correction is given by the Thouless-Valatin formula [2j| whose derivation in the context of the 1/N C - 
expansion is considered in Appendix iBl 

IV. TOY MODEL 

We want to study a toy quantum mechanical model imitating the SU(N c ) co \ OI ® SU(2) ap - m (g> SU(3)n symmetry of 
QCD. The Hamiltonian of this model is a generalization of the Lipkin-Meshkov-Glick model |27j : 



C fl J2,f[,f2 = l Sl,S2,si,S2 =1 \C = 1 
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xV f^' 2 f{s[\f 2 s 2 f lSl \J2^ 2S2C ^f lSl cj ■ (37) 

Here i/jf sc are fermionic operators ("quarks") with the "flavor" index / = 1,2,3, the "spin" index s = 1,2 and the 
"color" index c = 1,2,..., N c . They obey the standard anticommutation relations 

{V , / 1 s 1 ci,'i/'/ 2S2C2 } = S fx f 2 S SlS2 S ClC2 . (38) 

Introducing the compact multiindex notation 

A = (f,s); f = 1,2,3; a = 1,2, (39) 

we see that model <|37[) belongs to the class of models Q17p. 

The SU(2) sp i n ® SU(3)fi symmetry allows the following structure of V 

2 2 

where 

= ^.x^- a . ( 42 ) 



*fM\f a fi=WfLf*> (44) 
and 6 m „ are arbitrary coefficients. With this choice of we have 

#toy = fell (V'/Vc'V'/Vc') (i>f sc i>fsc) + &12 (^/V^/'v) (V'/.c^/'sc) 

+^21 (V/Vc^/'sc') {^Jso^fs'cj +h2 (V/Ve'V'/sc') f self's' c) ■ (45) 

V. EXACT SOLUTION OF THE TOY MODEL 

Let us diagonalize Hamiltonian 145fl . It is convenient to consider the pair of indices l|39|l as an SfJ (6) multiindex. 
As shown in Appendix IO the fermionic realization of the SU(6) Casimir operator is 



qSU{6) 

\ • ' us define the "fermion number" 



"spin" 



= ^V'/.c , (47) 



Ji = ^f s > c (^i)s's^fsc (48) 



and "flavor" generators 

1 



Ta = ^, sc {K)f'rtfsc (49) 
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where Oi and A a are Pauli and Gell-Mann matrices respectively. Using identity l|C.3(l from Appendix we find 



3 1 1 



(50) 



(51) 



Now the Hamiltonian (|4B|) can be rewritten in the form 



H 



toy 



1 

2~K 



+6 21 + +622 (2C. 



»=i 



-,S£/(6) 



6 



-iV. 



(52) 



This expression solves the problem of the diagonalization of the Hamiltonian. Indeed, all operators appearing here, 

8 3 



a 



SU(6) 



(53) 



0=1 



i=l 



commute with each other. 

We are interested in the "baryonic" color-singlet states containing N c quarks 



: cic 2 ...c Nc V'/ lSlCl V , ^ S2C2 ■ ■ ■ Tpfi 



+ 

iV c SiV c CN C 



|0>. 



For these states 



N 4 = N C . 



(54) 



(55) 



The states are totally symmetric in the SU(6) multiindices Ak — (fkSk) and correspond to the SU(6) irreducible 
representation of dimension 



D 



(N c + 5)! 
5!iV c ! 



with the eigenvalue of the Casimir operator 



C SU(6) = ^ Nc{Nc + 6) 



12 



(56) 



(57) 



This irreducible SU(6) representation has the following decomposition in terms of SU(2) sp - m ® 5J7(3)a irreducible 
representations : 



[Nc/2] 

>9 (fc) W(3) 



fc=0 



Here [iV c /2] stands for the integer part of N c /2. The SU{2) spins are given by 

jW = ^-fe (fc = 0,l,...,[7V c /2]) 

and the g^- 1 parameters of the 5C/(3)a irreducible representations are 

p (fe) = 7V C - 2k, q ik) = k . 



(58) 



(59) 



(60) 
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For the squared spin operator we obviously have 



(61) 



Using Eqs. 0, (EOJ), we obtain 

8 



a=l 



(p(*), 9 (*)) 



(62) 



Now we find the spectrum of energies of the toy model inserting Eqs. (|55|l . (|57(l . If) . I|62|l into Eq. (|52(l 



IiY c + l + AjW ( jW + i) 



AjW(jW + i ) + iiv c 



+ 622 (iV c + 5) 



(63) 



If we concentrate on the states with jW — 0(N®), then the large- N c structure of this exact result is as follows: 



E, 



(*) 



(fc) 



where 



E x = - 
1 2 



fen + ^22 + - (&12 + b 2 i) 



(64) 



(65) 



Eq = - (bu + 5622) 



(66) 



E ( _ k) 1 =(b 12 +b 21 )J^(J^ + 1). 



(67) 



The spins (|5§j) and the SU(3)n representations {jp^ k \q^) l|6"U|) of states, appearing in the spectrum of the toy 
model, coincide with the quantum numbers (|14|l of the nonexotic £ = states in the rigid-rotator model. Moreover, 

the dependence of the N~ 1 E < ^1 part of the energy (|67|l has the same J( J + 1) form as the spectrum of the rigid 
rotator (|15|l . However, the states with a nonzero "exoticness" (|1(J|) £ ^ are absent in the exact solution of the toy 
model. 

VI. LARGE- Af c APPROACH TO THE TOY MODEL: LEADING AND NEXT-TO-LEADING ORDERS 

In the previous section we extracted the l/iV c -expansion (|64|l for the spectrum of the toy model directly from the 
exact solution l|63|l . On the other hand, it is instructive to construct the systematic l/iV c -expansion for the toy model 
directly, without using the exact solution. The general formalism of the l/iV c -expansion for fermion systems was 
described in Section HTT1 Now we want to apply this formalism to the toy model l|45|l . 

We have to solve Hartree equations i|25|) . i|26|) . In principle, these equations have many solutions describing various 
branches of the spectrum of our model. We are interested in the solution corresponding to the baryonic states (color- 
singlet states with — N c ). In order to find this solution, let us make use of the analogy with the large- iV c QCD. 
In QCD the large- N c baryons are described by a mean-field solution (of some unknown effective theory) which has 
a combined SU(2) symmetry with respect to simultaneous space and flavor rotations. In the toy model we are also 
interested in the solutions of Eqs. I|25|l . I|26|) which are invariant under combined spin-flavor SU(2) rotations. Let us 
introduce projectors (k — 1,2,3) compatible with this spin-flavor symmetry: 



P 



/2S2/1S1 



f 2 6 SlS2 + Sf lS J Slf2 ) (1 - S fl3 )(l - S h3 ) , 



(68) 
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P hlf lSl = \ (*mJs iS2 - Sf lS2 5 Sl f 2 ) (1 - S fl3 )(l - 5 h3 ) , (69) 

Pi!lf lS1 =S SlS Jf l3 S h3 . (70) 
Using compact multiindex notation l|39|) . we can write 

J2p«Ml = s ki Pi% (7i) 

B 

j2 P AB=S AB . (72) 
fe=i 

The dimensions of these projectors can be read from their traces: 

SpP (1) = 3 , SpP (2) = 1 , SpP (3) = 2 . (73) 
Now let us consider the linear combination of these projectors 

3 



n AB = J2 £kP AB- (74) 



,(*) 

k=l 

Combining this with Eq. I|71|l . we see that ttab obeys the equation 

5> AB Pj^= £fe P^. (75) 

B 

Comparing this equation with Eq. I|25|) , we conclude that representation (|74|) is compatible with the Hartree equation 

We are interested in the baryonic states which correspond to occupying one level with N c quarks. The one- 
dimensional projector P^ is a perfect candidate for this occupied level. Therefore we take in Eq. (|24|l 

^ cc) -E^vr-^ 2 2- (76) 

i: occ 

Since the projector P^g is one-dimensional, the sum in Eq. ((75) contains only one occupied level: 

<t>?s = ^f* (£12 = -621 = 1). (77) 

Now equation leads us to the following result: 

kab = Vab\a<b>P { b} ai ■ (78) 

The leading order of the l/iV c -expansion for the energy E 1122(1 is given by the Hartree term (|23J) . Using Eq. (|76|) . we 
find 

E = ^P b \>Va>b>\abP b 2) a + 0(N° C ) . (79) 
With the explicit expressions for Va'B'\ab <|40(1 and for P^ (|69|) we easily find 

N, 



E : 

2 



b\l + &22 + g ( & 12 + b 2l) 



O{N c ) . (80) 



This coincides with the leading 0(N C ) part Ij65(l of the exact solution l|64|l . 

According to Eqs. (|71|l . (|74|) the single-particle energies are given by the following expression: 

_ Sp(ttPW) _ Pl^B^Pg^ . . 

fc s p p« SpP« ' { 1 
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Using the explicit expressions for (|68|) - l|70|) and for V (|40|l . we find 

£1=611 + 2(613 + 621), (82) 

£2 = 6n + 622 + ^ (612 + 621) , (83) 

£ 3 = 611 + ~6 2 i ■ (84) 

In the toy model (|45(l all RPA eigenfrcqucncies Q.^ vanish (see Appendix [DJ) so that the harmonic excitations 
are absent. According to Eq. Q36|l . in the absence of the harmonic excitations the 0(N®) correction Eq to the energy 
of the ground state is given by the equation 

E ° = -\ E E( £ ™- £ *)> ( 85 ) 

m:nonocc ?:occ 

where the summation runs over the single-particle occupied states £j and nonoccupied states e m with the weights 
corresponding to their degeneracy. In our case according to Eq. I|73[l the degeneracy of non-occupied levels £1 and £3 
is 3 and 2 respectively, whereas the occupied level £2 is not degenerate. Therefore we find from Eqs. lf%2")) . (f%3f . JBU, 

JEHJ 

E = -- [3(£i - £ 2 ) + 2(£ 3 - £ 2 )] = -622 + ^12 . (86) 
This result agrees with the expression for E (|66|l obtained from the exact solution. 



VII. ROTATIONAL EXCITATIONS IN THE l/iV c -EXPANSIONS 

As mentioned above, the mean-field solution ttab (|78[) violates both SU (2) S p ln and SU(3)fi symmetries but it is 
invariant under combined spin-flavor SU(2) rotations. Therefore the systematic 1 /Ac-expansion should involve a 
correct treatment of the collective coordinates corresponding to the rotation of the mean field. Since the spin rotation 
of the mean field is equivalent to the flavor rotations, it is sufficient to consider only the SU (3) flavor rotations. The 
invariance of the mean field under the As transformations leads to the J% constraint I0J . 

The dependence of the energy on the spin J appears in the order 0(N~ 1 ) and has the form 

A£., = ^!, (87) 

where / = 0(N C ). In the framework of the systematic l/A c -cxpansion the "moment of inertia" / can be computed 
using the large- N c version of the Thouless-Valatin formula [23( (see Appendix [Bj . Applying this formula to the toy 
model, we find [see Eqs. (|E.2|) . ljE.3|) in Appendix IE] 



/ _ V \ ' wi|A..|mi; \ T (0-1 



N c v-» ( (ii|A 3 |mi) \ , ^ / (m2|A 3 |i 2 ) 



2iZ2:occ mim2:nonocc 



TOl|A 3 |?l) J V Im 1 i 1 m 2 i2 I (i2|A 3 |m 2 ) 



Here S is the matrix (|34ll entering the large- N c RPA equation l|32|) . Strictly speaking, matrix S is degenerate so that 
expression (I88f) has to be rewritten in a more careful form. The details can be found in Appendix El The resulting 
expression for / in the toy model H45[) is given by Eq. i|E.6|) 



1 2(6 12 + 6 21 )- (89) 



Inserting this expression into Eq. 1)8 70. we obtain 



A^=^+^ J(J + 1) , (90 ) 



which agrees with the N c 1 part (|67|l of the exact expression for the energy. 
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VIII. PROBLEMS OF THE RIGID-ROTATOR APPROXIMATION 

The expression l|88|) for the moment of inertia I requires the knowledge of the RPA matrix S. The detailed 
description of this matrix is given in Section lTTTl [see Eq. (|35fl ]. If one omits the direct V contribution to S [represented 
by the W in the RHS of Eq. (|35|l ] then one arrives at the following approximation for S: 

Sniiii7ri2i2 ~ ^mim2^iii2 (^^l ^il ) ■ (9-Q 

We stress that this approximation is not justified by the large- N c limit. 

Inserting this approximate expression into Eq. (|88l) , one obtains the Inglis expression for the moment of inertia p4j 

I 1 = E^y2 ^ A3 l m ^ TO l A3 ^ (92) 



i:occ m:nonocc 



This expression can be also "derived" in the rigid-rotator approximation if one concentrates only on the rotational 
degrees of freedom instead of the systematic l/iV c -expansion (see Appendix IFl and Ref. (28|). In this way one arrives 
at the matrix i|F.4|) 

T N c ^ ^ (i\X a \m)(m\X b \i) + (i\X b \m) (m\X a \i) 

lab - ~t 1^ 1^ • (93) 

occ m: nonocc 

Here the occupied state i is associated with the projector P^ Q69[l. whereas the nonoccupied states m correspond to 
the projectors P^ (|68|l and P^ (|70|) . Due to the SU(2) spin- flavor symmetry of the mean field, the matrix I a i, is 
diagonal: 

h = hi = hi = hs , (94) 

h = hi = hb = I&6 = hi ■ (95) 

The transition £2 <-> £1 in Eq. 193|l contributes to the moment of inertia h, whereas the transition £2 «-> £3 generates 
h- We find from Eq. 



N c ^ S P (P«A a P( m )A h + P«A fc P(" l )A Q ) 
^ = T" 1. 7— T- 



4 

m— 1,3 



(96) 



which results in 



h - 2(-2fe - M ' (9S) 



In the rigid-rotator approximation the spectrum of low- lying 0(N C x ) excitations is described by Eq. I|15|) : 

J(J+1) 
2/i 



(99) 



This expression has the same structure as Eq. (|87|) derived in the systematic l/A^-expansion. However, comparing the 
rigid-rotator moment of inertia h H97|) with the correct expression / (|89|l . we see that the rigid-rotator approximation 
gives a wrong result: 

h+h (100) 

The reason is obvious: approximation l|91l) is not justified by the large-Ac limit. Approximation l|91ll has the same 
problem in the case of the model with the flavor group SU(2), see Appendix IU1 Since the Inglis formula l|92l) fails 
already in the nonexotic sector of 0(N~ 1 ) excitations, there is no sense to discuss the results based on the Inglis 
expression for h- Moreover, as mentioned above, the states of the rigid-rotator Hamiltonian Q with a nonzero 
"exoticness" ()l(Jfl £ ^ are simply absent in the spectrum of the exact solution of the toy model l|45|) . 
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IX. LARGE-TV, MODELS WITH COLOR-SINGLET BOSON FIELDS 

The large-iVc counting in QCD is compatible with the large- N c behavior in various models: naive quark model, 
Skyrme model 0,0], Nambu-Jona-Lasinio model (NJL) [23, chiral quark-soliton model (CQSM) j2g etc. Some 
of these models (e.g. NJL) contain only colored fermion fields, others (e.g. Skyrme model) involve only color-singlet 
bosonic fields, some models such as CQSM contain both color-singlet boson fields and colored quark fields. Let us 
concentrate on two extreme cases of models with local (in time) actions containing 

1) colored fermion fields, 

2) color-singlet boson fields. 

In spite of the difference between these two groups, the analysis of the l/i\T c -expansion has many common features. 
Indeed, we can bosonize the interaction between the fermions and integrate out the fermions (see Appendix EJ. In 
this way one arrives at a theory with color-singlet boson field it but the action S[n] of this theory is nonlocal. 

Whatever action S[ir] we have, local or nonlocal, the spectrum of 0(N®) excitations 

E = const + ^n v n v + 0{N- r ) , n v = 0, 1, 2, . . . (101) 

V 

is given by frequencies u> = f2„, which are determined by the equation 

detJfH = 0. (102) 
Here K(u>) is the quadratic form of the boson action S[tt] in the frequency representation: 

K{u>) = ±r [ dte-*" t / / '?L ■ (103) 

V J N C J STT(t)6TT(0) V ' 

In theories with nonlocal actions S[n], obtained via a bosonization of local theories with colored fermions, one can 
reduce equation (|102(l to RPA equations as shown in Appendix 1X1 

In theories with local actions S[ir] quadratic in time derivatives, we have the following structure of the quadratic 
form K(u>): 

Kluj) = -K 2 u? + iK x uj + K . (104) 

In this case Eq. (|102|) is equivalent to the existence of a periodic solution of the equation 

(K 2 d 2 t + Kidt + Kq) Sn{t) = (105) 

with the frequency u>. This equation has a simple physical meaning: it describes periodic solutions of the classical 
equations of motion corresponding to an infinitesimal deviation from the static solution. In a general case this problem 
cannot be reduced to the diagonalization of the rigid-rotator Hamiltonian. 

However, in Ref. |2(J (see also Ref. |2lJ) a simple model was considered whose spectrum of 0(N®) excitations 
contains a branch which can be successfully described by the rigid-rotator approximation. This model involves a 
particle with mass fi and electric charge e, moving in the 3-dimensional space in the field H of a monopolc with 
magnetic charge g 

H = g-r^ = V x A + [string field] . (106) 
l x l 

In addition the model has a central potential V so that the total Lagrangian is 

L=^!-e(xA)-nW). (107) 

If one interprets this Lagrangian as a sort of effective theory for large- N c "baryons" coming from some "microscopic" 
theory, then the large- N c limit corresponds to the semiclassical limit for the Lagrangian l|107l) : 

fx = 0(N c ), V = 0(N c ), (108) 
eg=^N c B, B = 0(N°). (109) 
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Here the integer number B has the meaning of the "baryon charge" in the underlying "microscopic" large- N c theory. 
The representation (|109() for eg agrees with the Dirac quantization rule for the magnetic charge. 

As shown in Ref. j2£|, this large- N c inspired semiclassical limit justifies the reduction of the model to the rigid- 
rotator Lagrangian 



dt 



egSl 3 



(110) 



for a subset of 0(N®) excitations. Here fi a is the angular velocity. The moment of inertia I is given by the equation 

/ = [j,R 2 , V'(R) = 0. (Ill) 
The Lagrangian l|110fl corresponds to the Hamiltonian 

1 



Hrot = ^j( J i + J z) 



(112) 



with the quantization constraint 



J'3 = eg = 



N r B 



(113) 



Here Jg is the projection of the momentum J on the third axis in the "body-fixed" frame. The spectrum of J is 



j = Nc\B[ N e \B\ 1 N C \B\ | 



As a result the spectrum of energies is 



E 



21 



j{j + i)-(4Y 



i 

21 



J(J+ 1) 



N r B 



By analogy with Eq. (|10|) one can introduce the "exoticness" 



N C \B\ 

£, = J 7; — = 0, 1, 2, 3, 



Then one finds from Eq. (|115fl 



E = Ee = — 
5 21 



NAB\ 



(2£ + !)+£(£ + !) 



In the limit I ~ N c — > oo we have equidistant 0(N®) excitations parametrized by £: 



4/ 



-(2£ + l) + 0(l/Ay. 



(114) 



(115) 



(116) 



(117) 



(118) 



As shown in Ref. [23] , the "exotic" spectrum (jll8[) of the rigid rotator l|110(l coincides (up to 1 /N c corrections) with 
a part of the true spectrum of the original model (|107|) . This observation was used as an argument in favor of the 
validity of the rigid-rotator approach to exotic states. 

Let us check that the spectrum Ijll8|l is described by one of the frequencies lu determined by the general equa- 
tion 1|102|) . The full spectrum of the harmonic excitations can be found from the solutions of equation Ijl05|l . Due 
to the symmetry of the problem, the equation of motion l|105|) can be split into two parts. In the case of radial 
infinitesimal oscillations Sx ra d(t) controlled by the central potential, one can ignore the K\ term associated with the 
magnetic field: 



(K 2 d? + K ) Sx rad (t)=0. 



(119) 



We also have the periodic motion Sx C i rc (t) along an infinitesimal circle due to the magnetic force. In this case we can 
omit the term K containing the central potential: 



{K 2 d 2 t +K 1 d t )5x clrc {t) = Q. 



(120) 
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In the frequency representation Eq. (|120f> takes the form 

( k ~lj 2 K 2 + iujK x ) Si clrc = . (121) 

This equation has nonzero solutions 8x C i rc for two frequencies to. The zero frequency to = corresponds to the J3 
degeneracy of the quantum spectrum. The second nonzero frequency 

u = — (122) 
A* 

describes the circular motion of the particle with mass /i and charge e in the magnetic field H. This frequency is 
responsible for the "exotic excitations" (|118(1 . Indeed, the magnetic field of the monopole is H = g/R 2 , and the 
moment of inertia I is given by Eq. I|lll|) so that 



Using Eqs. ifTMjl . fTO . (tTTTIt . we find 



w=y- (123) 



- = ™ = O(iV c ). (124) 
This frequency agrees with the spectrum of harmonic excitations H118(l : 

E i+1 -E i = ^ + 0(l/N c ). (125) 

In addition one has excitations determined by the eigenfrequency of Eq. I|119|) ■ 

Thus in the case of the simple model ()107(l the rigid-rotator approximation leads to the spectrum which is a 
part of the full 0(N®) spectrum determined by Eq. i|l(J2|) . However, this agreement of the rigid-rotator spectrum 
with the general equation (|1(J2|) is "accidental" . In a general case one has to solve equation (|102f) or the equivalent 
problem (|105fl in order to find the spectrum of 0(N®) excitations and this solution has nothing to do with the rigid 
rotator approximation. In Appendix [H] we describe an extension of the model p()7|l which explicitly demonstrates 
the breakdown of the rigid-rotator approximation. 

Let us summarize. Both in the models with colored fermions and in the models with color-singlet bosons the 
spectrum of the 0(N®) excitations is determined by the solutions of equation (|102fl . In a general case the rigid- 
rotator approximation cannot describe these excitations. 



X. CONCLUSIONS 



In this paper we have considered several simple quantum-mechanical models imitating the large- N c QCD. The main 
aim was to understand the status of exotic baryons at large N c . Our results confirm the conservative opinion that 
the exotic states are only a part of the 0(N®) excitations whose properties are determined by the quadratic form of 
the action written in terms of color-singlet boson fields. Recent attempts to extract reliable information about exotic 
baryons directly from the rigid-rotator Hamiltonian are incompatible with the results of our analysis. 

We have considered two different classes of models. The first group includes models containing colored fermions. 
The second class of models is based on local actions for color-singlet fields. In both cases the spectrum of the 0(N®) 
excitations is determined by the eigenfrequencies of the quadratic form of the action. But in the case of colored 
fermions one has to work with the nonlocal bosonized action. Our analysis of this nonlocal quadratic form was based 
on the reduction of the problem to RPA equations. In the case of models with local actions for color-singlet bosons 
it is convenient to perform the analysis of the quadratic form in terms of the semiclassical quantization of periodic 
classical solutions. 

We have studied an exactly solvable toy model with colored fermions. We have shown that in the large- N c limit 
this model is described by a mean-field solution whose spin-flavor symmetry coincides with the symmetry of the 
large- N c baryons in QCD. Using this mean-field solution as a starting point for the systematic l/A c -expansion, we 
have computed the energy spectrum in the first three orders of the l/A c -expansion. The results of the 1/ A c -expansion 
are in complete agreement with the exact spectrum of the model. 

We have also tested the rigid-rotator approximation. The results coming from the rigid-rotator approximation are 
rather disappointing. The effective rotational Hamiltonian derived in the context of the rigid-rotator approximation 
has exotic baryonic states which are absent in the exact solution of the model. 
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The rigid-rotator approximation also has problems in the nonexotic sector. Although the rigid-rotator approxima- 
tion predicts the correct quantum numbers for the nonexotic low-lying 0(N~ ) baryonic excitations, the values of the 
energy coming from the rigid-rotator approximation are wrong. The comparison of the rigid-rotator approximation 
with the systematic l/iV c -expansion allows us to detect the origin of the problem: the rigid-rotator approximation 
leads to the Inglis formula which ignores the quadratic form of fluctuations around the mean field. On the contrary, 
the calculation of the rotational spectrum based on the systematic l/A c -expansion results in the Thouless-Valatin 
formula which involves the coefficients of the RPA equations. 

Although our model with colored fermions reproduces the main symmetries of QCD, at one point there is an 
important difference: the absence of the exotic states in the toy model is simply a model artifact: the "baryonic" 
states in the model contain only N c quarks above the bare vacuum. This leaves no possibility to have extra quark- 
antiquark pairs which could generate exotic quantum numbers. In QCD the situation is different: the existence of 
exotic states is allowed. However, the question is whether the rigid-rotator approximation can correctly describe the 
exotic baryon states in QCD. Our analysis of the toy model raises doubts about the relevance of the rigid-rotator 
approximation. Actually it is not a problem to construct solvable models containing exotic states and to show explicitly 
that the rigid-rotator approximation does not describe properly these exotic states. 

We also have considered a group of models with local actions for color-singlet boson fields. This group contains a 
simple model compatible with the rigid-rotator approximation. However, our analysis shows that this agreement is a 
consequence of the simplicity of the model. A slight extension of this model leads to the breakdown of the rigid- rotator 
approximation in the exotic sector. 

To summarize, our results agree with the traditional approach to the 0(N®) baryonic excitations (including the 
exotic ones) as states controlled by the quadratic form of fluctuations of the effective action. Within the framework 
of the large-7V c expansion in QCD (without extra model assumptions) we see no justification of the rigid-rotator 
approximation. 

Certainly, one should be rather careful with the extension of our toy-model experience to the case of the large- N c 
QCD. For example, the fermionic toy model is based on the quantum mechanics of spin systems in O-dimensional 
space. This model has colored fermions but no gluons and no gauge symmetry. In the model calculations we were 
interested in the discrete eigenstates of the Hamiltonian, whereas in reality exotic baryons are resonances. One also 
should keep in mind that this toy model imitates the strict chiral limit of QCD. In the real world the situation is 
more difficult since one has to deal with a nontrivial interplay between two small parameters: the light-quark mass 
matrix and 1/N C . 

Acknowledgments. The author appreciates many interesting and instructive discussions with Ya. I. Azimov, 
T. D. Cohen, D. I. Diakonov, V. Yu. Petrov and M. V. Polyakov. This work was supported by DFG and BMBF. 



APPENDIX A: PATH INTEGRAL DERIVATION OF HARTREE AND RPA EQUATIONS IN THE 

LARGE- iV c EXPANSION 



In this appendix we use the path integral approach in order to construct the l/-/V c -expansion for models with 
four-fermionic interaction and show how Hartree and RPA equations appear in the context of the l/A c -expansion. 
We start from the Lagrangian corresponding to the model Hamiltonian l|17fl 



L = itfj\ c d t 1pAc ~ {^A'c'^B'd) VA'B'\AB (lp~A c lpBc) 



(A.l) 



One can bosonize the four-fermion interaction, introducing auxiliary variables ttab and integrating out fermion vari- 
ables ip, ip + : 



Dtp / Dip + expii / dt 



itpi c d t 1pAc ~ Oa'c'V'B'c') V A 'B'\AB (^Ac^Bc) 



Dtt / Dip / Dtp + expli / dt 



ilpA c d t lpAc - IpAc^ABlpBc + -y7rA'B'(V l ) A' B' \AB^ AB 



Here we have the effective action 

S[ir] = N c 



= J Direx~piS{iT} . 
-ilndct(i<9 t - n) + - J dtv A > b\V~ X ) A' b>\ab-k ab 



(A.2) 



(A.3) 
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The presence of the factor N c in this action allows us to apply the saddle point method for the construction of the 
l/-/V c -expansion. The saddle point equation is 

' !>S;!J =i{t,B\—^\t : A) + (V- 1 ) AB \ A>B> n A , B , =0. (A.4) 



N c 5lT AB (t) ' id t -TT 

For static saddle points it we have 

where <fig and e n are eigenvectors and eigenvalues of the matrix n : 

= e n ^ } . (A.6) 
The ±i0 pole uncertainties in Eq. (|A.5(I should be resolved relying on the physical picture of occupied states 



£E^- = E^-^ C) . (A.7) 



After the oj integration, the summation runs only over the occupied single-particle states (p^K 

Inserting Eqs. (|A. 5|) and (|A.T|) into the saddle-point equation (|A.4(1 . we obtain the Hartree equation 

kab = V ABlA , B ,P B °, A ) . (A.8) 
Computing the action (|A.3|I for this ir AB , one easily finds the 0(N C ) Hartree energy 

E = ^V m A'B>P$?P$$ + 0(N° C ) . (A.9) 

Now we want to compute the 0(N®) corrections. In a general case these corrections will describe the harmonic 
excitations with some frequencies f2„ 

AE = ^nM» = 0(N°), (A.10) 

where numbers n v = 0, 1, 2, . . . are integer. In order to find these frequencies we must consider the quadratic form of 
the action (|A.3|I in the frequency representation: 

""'"^/'■'"- ......nLw ^ 1 "' (A ' n) 

The spectrum of the harmonic excitations (|A.10fl is determined by the frequencies u> = tt v at which the quadratic 
form K becomes degenerate: 

detK(Q v ) = 0. (A.12) 
The calculation of the quadratic form of the action (|A. 3|) is straightforward: 

1 6 2 



N c 8nA 1 B 1 (ti)5nA 2 B 3 (h) 
I „ , 1 



S[tt] 



= t(h,B 2 \--^—\t u A 1 )(t u B 1 \--^—\t 2 ,A 2 ) + (V- 1 ) A2B2 6(t 1 -t 2 ). (A.13) 
Now the quadratic form (|A.11|I can be represented as follows: 

Ka 1 b 1 \a 2 b 2 (u) = iL AlBl \A 2 B 2 (^) + i v ^) AlBl \ A2B2 > ( A - 14 ) 



where 



. ULii , ;S » = / dtie-** /(0,P 2 |^— -\t u A 1 ){t 1 ,B 1 \^— -|0,A 2 ) 



*E E 

z:occ m:nonocc 



(fl2|i><*|Ai)<fli|m)<m|A 2 ) , (B a |m)H^i><Si|i)<i|i4 a 



e m - £i - w - »0 



£ m - £i + a; - zO 



Here we use the standard Dirac notation 

(A\k) = (k\A)* = . 
Taking some antiunitary operator (e.g. time reversal) 



and introducing compact notation 



7r e = 07TIT 1 



|AB) - \A)®\B). 



D = 7r®l-l®7r, 



we can write 



(il^XBilm) = (fli|m)(A?|i fl ) = (Pu^rm*) 



Then expression ljA.15|) takes the form 

^A 1 B 1 |A 2 B 2 (w) = 

Here we have introduced projectors 



D\mi e ) = (e m -e l )\mi 9 ) 



Pi = J2 E \mi 9 ){mi 6 



t:occ mmonocc 



p 2 = j2 E i imfl ><* 



Next we define 



t-.occ m:nonocc 



PtP - ShPi 



E = P x - P 2 , 



Q(w) = E- 



£> - w - iO • E 
Then the RHS of Eq. (|A.23|) takes the simple form 
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Inserting this into Eq. I|A.14|) . we see that 



l A 1 B 1 \A 2 B2 



(A.30) 



Although 9 is an antiunitary operator, for any given basis {|^4)} we can find such a unitary operator U that (for the 
basis states \A) only) 



\A e ) = U\A) . 

Now expression (jA.30|) becomes 

K AxBx]A2B M = (B 1 A 1 \(1 ® f/ + )QH(l ® U)\A 2 B 2 ) + {V- l ) AiBM2B2 
Next we define the permutation operator II 

H\AB) = \BA) . 

This allows us to rewrite relation i|A.32|) in the operator form without indices A^ , B^ 

K{u) = n(l ® U + )Q(oj){\ ®u) + v- 1 . 

Let us define 

K(u) = (1 ® J7)nJC(o;)(l ® [/+) , 



(A.31) 
(A.32) 
(A.33) 
(A.34) 
(A.35) 



Now it follows from Eq. (|A.34|) 



The next step is to define the projector 



V= (1 ® U)VU(1 ® Z7+) . 



P = Pi+P 2 , 



(A.36) 
(A.37) 
(A.38) 



and to notice that 



P z = P 



= PQ = QP = PQP. 



(A.39) 



(A.40) 



The operator Q can be inverted on the subspace of the projector P. Let us consider operator P 1 + V ) P on 11 lis 
subspace and let R be its inverse: 



R 



P 



(q- 1 



V) P 



PiQ^ + VjP 



R = P, 



(A.41) 



RP = PR = R. 



For brevity we will write 



R = 



P(Q- 1 + V)P 



= P 



PVP) P. 



(A.42) 



(A.43) 



implying that all ambiguities should be understood in the sense of Eqs. I|A.41(I . I|A.42J) 
Multiplying equality ljA.41|) by Q and using Eqs. l|A.40fl . (|A.42|i . we find 



R + QVR = Q . 



(AAA) 
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Using Eq. (|A.37|) . we can write 

k (v - vrvJ = {q + v- 1 ^ (v - VRV^j =i-rv+(q- qvrJ V . 

We can simplify the RHS using Eq. 1A.44I) : 

K ( V - VRV 1=1. 



(v - VRV^j 



Combining this with Eq. (|A.43|) . we find 

R- 1 = V - VRV = V - VP (Q- 1 + PVP ) P\ . 
Inserting Eq. I|A.28|) . we arrive at 



£L> + PVP - wS - iO 



= v -vp 

According to Eq. i|A.35|l we have 

\det K(lu) \ = \det K(lu)\. 
Therefore Eq. l|A.12fl is equivalent to the equation 

detK(Q u ) = 0. 

Combining this with Eq. I|A.48|I . we see that the problem reduces to the equation 

dct (ZD + PVP - = , 



PV . 



(A.45) 



(A.46) 



(A.47) 



(A.48) 



(A.49) 



(A.50) 



(A.51) 



where the determinant must be computed on the subspace of the projector P. This is equivalent to the eigenvalue 
equation 



SZ V = n v T,z v 



where 



S = £L> + PVP . 



(A.52) 



(A.53) 



Equation (|A.52j> is nothing else but RPA equation lO- Indeed, according to Eqs. (|A.20|) . HA. 221) . 1A.24|) . (|A.25jl . 
(|A.27jl we have 

T,D — /Z {\mi e }(mi e \-\im e )(im e \) (tt (g) 1-1® n e ) 



f.occ m:nonocc 



Yl (£m-£i)(\mi e ){mi e \ 

i:occ m:nonocc 



9 ><W|) 



(A.54) 



Using Eqs. (|A~24)l . (|A~25)l . (JOTJ), l|A~36|) . l|A~38)) . we can write 

PVP= J2 J2 {\mi e )(mi e \ + \im e )(im e \) 



z,j:occ m,n:nonocc 



X(l ® U)VH(1 <g> C/+) + bn e )(jn e |) 



J] 5^ (|rm 9 )(rm| + |im 9 )(H) V (|jn)(n/| + |nj)(jn 9 |) . 



i,j:occ m,n:nonocc 



(A.55) 
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Now we use Eq. <|31[) 

PVP= ]T E [\mi°)W mijn (nj \ + \im B )W imjn (nj e \ 

i,j:occ m,n:nonocc 

+ \mi 9 )W mmj (jn e \ + \im 9 )W lmnj (jn 9 \] . (A.56) 

Taking the vector Z v in the form 

Z U = J2 E {X mt \mi e ) +Y rnl \im 9 )) , (A.57) 

i:occ m:nonocc 

we find using Eqs. lA.27jl . (|A.54jl . (|A.56|) 

sz ' y = E E (a^O-^™ 9 )) , (a.58) 

z:occ m:nonocc 

(S£>)^ = ^ ^ ( £m - £4 )(^l™ 9 )+^l^» , (A.59) 

i:occ mmonocc 

(pUP)iT = E D™') (^ J ni + W minj Y^) 

2,j:occ m,n:nonocc 

+ |zm e ) (W imin A%- + W imnj Y^)] . (A.60) 
Inserting these expressions into Eqs. (|A.52I) . 1A.53|) we obtain 

(e m -e i )X^ i +Y J E { W ™jnX v nj + W minj Yr)=n v X» mi , (A.61) 

j:occ n:nonocc 

(e m -s l )Y^+J2 E {W imj nXn j +W imnj Y» j )=-n v Y» i . (A.62) 

j:occ nmonocc 

Using relations 120|) . I|21|) . we see that these equations coincide with RPA equations (|27H . (|28l) . 

APPENDIX B: THOULESS-VALATIN FORMULA AT LARGE N c 

In this appendix we derive the large- N c version of the Thouless-Valatin formula |23| for the moment of inertia 
controlling 0(N^r 1 ) rotational excitations (|87J) . This formula is well known in the many-body physics 22] . The 
derivation presented below follows the traditional ideas but emphasizes the role of the large N c in the justification 
of Thouless-Valatin formula. We start from a rotationally invariant Hamiltonian H commuting with the angular 
momentum Jk 

[H,J k }=0. (B.l) 

In the case of the rotational excitations l|87l) the dependence on the spin J appears in the order 0(N~ 1 ). In principle, 
the calculation of the NNLO 0(N C T 1 ) corrections to the energy requires a serious work. But if one is interested only in 
the J dependent part of this 0(N~ 1 ) contribution, then the result can be computed using Thouless-Valatin formula. 

Let us consider the auxiliary Hamiltonian H — ujJ% instead of the original Hamiltonian H and let us study the 
ground state energy E^^uj) of the Hamiltonian H — luJ^ as a function of u in the interval 

N~ l < uj < iV° . (B.2) 

In this intermediate region we can treat ui either as a "numerically small" 0(N®) quantity or as a "numerically large" 
0(N~ 1 ) quantity. Both approaches allow us to derive corresponding expressions for the ground state energy E^°\w). 
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Matching the results in the intermediate region (|B.2|I , one can obtain a simple representation for the moment of inertia 
I controlling the J( J + 1) splitting (|S7J> . 

If we consider the case of ui = 0(N®), then the ground state energy E^°\uj) of the Hamiltonian H — uiJ^ is given 
by the u> dependent Hartree energy iVc-Ei (u>) 



= NeEi(u) + 0(N°) . 

Let us compute the second order derivative of this expression with respect to u> in the interval ljB.2 



duo- 



UJ 



d 2 



N B Ei(w) 



d 2 



duo' 



rJVc-Bl(w') 



(B.3) 



(B.4) 



'=0 



Now let us change the approach and consider w as a "numerically large" 0(N C : ) quantity. We start from the 
l/iV c -expansion for the Hamiltonian H (without the C0J3 correction) 



Ef = N c Ex +E a + — [ Cl J(J + 1) + 02] + 0(N~ 2 ) . 



(B.5) 



According to Eq. IB. 1(1 Hamiltonian H commutes with J3. Therefore we find from Eq. I|B.5|) for the energy E^°\u)) 
of the ground state of the Hamiltonian H — ujJ^ 



E^(u) = minf 
J, J3 V 



min I N c Ei +E + — [d J(J + 1) + c 2 ] - w J 3 + 0(W C ~ 

J,J3 I iV r 



Obviously the minimum is achieved at J3 = J so that 



E {0) (lu) = min \n c Ei + E + -L [ Cl J(J + 1) + c 2 ] - w J + 0(7V,T 2 ) 
J I iV c 



The point of the minimum is 



J = 



uiN c — ci 
2^ 



(B.6) 



(B.7) 



(B. 



Due to condition (IB.2II we have 



Now we find from Eq. (|B.7I 



1«J*^«N C . 



(B.9) 



We see that 



£ (0) (cj) w NcEi+Eo 



4 Cl 



9w 2 



AT, 



Combining this result with Eq. I|B.4|) . we see that 

<9 2 



&y 2 



AT c £i(u/) 



w'=0 2C 1 



(B.10) 



(B.ll) 



(B.12) 



Comparing Eqs. (|B.5|) and H87fl we find the connection between the coefficient C\ and the moment of inertia J: 



ci 



No 
21 



(B.13) 
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Now we find from Eq. (|B.12|I 



1 = ~ N - ^ 



(B.14) 



This intermediate result simplifies the calculation of the moment of inertia /. Indeed, in the beginning the problem 
of the calculation of / was a NNLO problem. But now it is reduced to the expression l|B.14() containing only the Hartree 
energy N c Ei (10) for the Hamiltonian H — lo J3, i.e. to the problem of the leading order 0{N C ). The calculation of the 
second order derivative with respect to lo in Eq. I|B.14I) is equivalent to the construction of the second perturbation 
theory in lo for the Hartree equation. 

Now we have to write the Hartree equation for the Hamiltonian H — lo J3 with 



^3 = ^ + ^3ip , 



(B.15) 



where 03 is the Pauli matrix acting on the spin indices. Due to the lo J3 term the expression for the Hartree energy 1231) 
is modified as follows 



£1 = \v. 



AB\A'B'PbA 



(occ) p(occ) ^ 



p(occ) 
AB BA ■ 



(B.16) 



Let us differentiate this expression with respect to lo: 



1 



dE 1 = hAb dP B T ] 

dtu duj 2 



(occ) 



{v^AB P BA 



(B.17) 



Here 



h A B - V A b\A'B>Pb'A' ~ 77 ( ct 3)ab 



(B.18) 



is nothing else but the one-particle Hartree Hamiltonian with the single particle eigenfunctions 



h-AB^ 



(B.19) 



and with the projector onto the occupied states 



P, 



(occ) 



BA 



E *W 

i: occ 



(B.20) 



The projector P B ° A depends on the parameter lo of the Hamiltonian H — L0J3. The standard perturbation theory 
allows us to compute the derivative 



dP, 



(occ) 



BA 



duo 



E 



E 



i:occ m:nonocc 



St - Sr. 



dio 



{m\ — \i)4> B 4>\ + {i\ — \m^ L; , 



(B.21) 



Combining this with Eq. I|B.19|) . we see that 



hAj 



dP, 



(occ) 



BA 



= 0. 



(B.22) 



We insert this result into Eq. (|B.17|I 



dE x 1 



dio 



~2 ^ a ^AB Pba ] - 



(B.23) 



Next we take the second lo derivative 



d 2 E 1 1 



dP, 



(occ) 



9a; 2 " 2^ )ab 8u 



(B.24) 
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Using Eq. (|B.21(1 . we find 



d 2 E 1 



2:occ m:nonocc 



dh 

{rn\ — \i){i\(T3\m) + (i\ — \m)(m\a 3 \i) 



According to Eq. (|B.18|) we have 



o p(occ) 1 
ur B'A' 1 



Inserting expression 1|B.21|) . we find 

dh AB 1 Y^ Y" 1 



9u> 



j:occ n:nonocc 



Now we take the matrix element of this equation 



E E r4 



0ft 



WmiinH^jli) + W minj (j\ — \n) 



dui 



j:occ n:nonocc ■> 

In the case when the state i is occupied and m is nonoccupied, let us introduce the notation 



2_ dh 



2 dh 



Then we find from Eq. (|B.2 



(e m - £j) a ml = (m\a 3 \i) — ^ ^ [W mijn a n j + W m injb n3 

j:occ n:nonocc 



Exchanging i «-> m in Eq. (|B.28(1 . we obtain the second equation 

(em - £t) = (i|o- 3 |m) - (w imjn a nj + W lmnj b n 

j:occ n:nonocc 

Obviously these equations can be rewritten in the form 



EE*. 

i' :occ 

where S m i m 'i' is the RPA matrix ()35D . 



z :occ m .nonocc 



(z|cr 3 |m) 



Now we find from Eq. I|B.25|) 



d 2 E 1 
duo 2 



]E E ( 

i:occ minonocc 



a m i(«|cr 3 |m) + b mi (m\a 3 \i) 



Combining this result with Eq. 1B.14|) , we obtain the moment of inertia 



I = ^ ^ (a m j(i|cr 3 |m) + b mi (m|o- 3 |i) 



2:occ m:nonocc 



(B.25) 



(B.26) 



(B.27) 



(B.28) 

(B.29) 
(B.30) 

(B.31) 
(B.32) 

(B.33) 



The result l|B.35|) for the moment of inertia / coincides with the well-known Thouless-Valatin formula |2j 
the subtleties which differ our large- N c matrix S (|35|l from the standard RPA matrix. 



(B.34) 

(B.35) 
I up to 
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APPENDIX C: 517(6) CASIMIR OPERATOR 

In this appendix we derive expression (|46() for the SU(6) Casimir operator. The N 2 — 1 generators t a of SU(N) 
can be normalized by the conditions 

S P (t a t b ) = ±S ab , (C.l) 
S P t Q = 0, (C.2) 

E ^cd = \ (sadSbc - ^S AB S CD ) . (C.3) 
o=l ~ x 



l AB l CD — n 

In the case of SU (6) the last relation results in 



1 /. . 1 



2 V 6 

Therefore the fermionic realization of the SU (6) Casimir operator is 

35 6 No 



E = ^ 5 AD S BC - -SabScd ■ (C.4) 



C 2 U{6) = E E E ^tc'tAB^B^mjiDM 



a=l A,B,C,D=1 c,c'=l 



" ^ ! I 

2(^Xc'^B,c')(^S,c^,c) - J^(V , Xc'V , A,c')(V , B,cV'B > c) 



E E 



(C.5) 



A,B=1 c,c'=l 

Returning to the (/, s) form i|3l3|) of the SU(6) multiindices, we arrive at Eq. (JIB}. 

APPENDIX D: RPA EQUATION IN THE TOY MODEL 

In this appendix we describe the structure of the RPA equation for the baryonic states of the toy model H45fl . In 
the Hartree picture for these states we have one occupied nondegenerate level <|83|) and two nonoccupied levels H82[) . 
(|84|l with degeneracies 3 and 2 respectively. This structure of the levels can be easily understood in terms of the 
unbroken spin-flavor SU (2) symmetry. The generators of this symmetry are 

G a = J a + T a (a = 1,2, 3) (D.l) 

where J a is spin and T a is isospin. The occupied state has G = 0, whereas the nonoccupied levels have G — 1/2 and 
G = 1 with the corresponding degeneracies 2G+ 1. The RPA equation also can be split into two sectors corresponding 
to G = 1/2 and G = 1. The nontrivial part of the RPA equation corresponds to the G = 1 sector. The RPA matrix S 
can be computed using the general expressions l|3~5l) . (|3"T1) . the explicit form of Va 1 a 2 |a 3 a 4 m the toy model (|4*U|) and 
the results for the single-particle levels (|82|l . (|83|l . In the G — 1 sector the RPA matrix S looks as follows 

Smini = 7: (612 + b 2 i) S mn (11) (m, n : G = 1) . (D.2) 



The precise form of this matrix depends on the phase conventions for the single-particle states. We fix the phases 
using the condition 

(m\X a \t)\_ x (I 
(i\Xa\m) 



= S am 1 J (a,m = l,2,3). (D.3) 



Here we label the states \m) of the G — 1 nonoccupied level with numbers m — 1,2,3. 

Inserting the result for the matrix S l|D.2(l into the general RPA equation l|32|l . we see that the RPA equation has 
only zero eigenfrequencies = in agreement with the absence of the harmonic excitations in the exact spectrum 
of the toy model. 
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APPENDIX E: ROTATIONAL EXCITATIONS IN THE TOY MODEL 

In this appendix we compute the moment of inertia which determines the spectrum of rotational excitations of 
baryonic states in the toy model 145|) . using Eqs. 1B.33|) . i|B.35(l derived in Appendix IB! 

In models, where the solution of Hartree equations has the spin-flavor symmetry, the projector p( occ ) onto the 
occupied states commutes with the generators of the combined spin-flavor rotations: 

[P (occ) ,K + r a )] =0. (E.l) 

Here the Pauli matrices a a , T a act on spin and flavor indices respectively. This symmetry allows us to obtain an 
equivalent representation for the moment of inertia by replacing 03 — > — T3 in Eqs. l|B.33jl . I|B.35(I . In the case of 
the SU(3) flavor group we must use the Gell-Mann matrix A3 instead of T3. As a result we arrive at the following 
representation for the moment of inertia 

/ = If ^ (a m i(i\*3\m) + b mi {m\\ s \i)) , (E.2) 

i:occ m:nonocc 

where a m i & n d b m i are solutions of the equation 

V" c ( a m'i> \ - f ( m \h\i) \ ( V o\ 

2^ b ™™'* I h m , v J - I (i|A 3 |m) I ' ^ 

i :occ m .nonocc ' 

and S is the RPA matrix 

We can rewrite Eq. (|E.2p in the form l|88|l but one must keep in mind that matrix S is degenerate. Using 
expression (|D.2I) for S m i m 'i' in the toy model, we find from Eq. I|E.3|I 



rn' — l 



2 

The solution of this equation is 

a m i \ = 1 ^ ( 1 + c 

bmi J bi2 + &21 V 1 — C 

where c is an uncertain constant. This constant cancels when we insert Eq. l|E.5Jl into Eq. (|E.2|I : 



(E.5) 



I = ~y y j (a ml (i\\3\m) + & TO i(m|A 3 |i)) = — -(a mi + b mi ) 



4 

m— 1 



Nr. 



2 b 12 +b 2 i 



(E.6) 



APPENDIX F: RIGID-ROTATOR LAGRANGIAN 



In this appendix we describe the structure of the effective rigid-rotator Lagrangian for models containing colored 
fermions. Let ttq be a static solution of the Hartree equations i|A.6|) . (|A.8(1 . If R is a transformation from the symmetry 
group G of the theory, then RitqRT 1 is also a solution of Hartree equations. Now let us apply a time-dependent 
transformation R(t) to the static solution ir Q 

tt^^R^ttoR- 1 ^). (F.I) 

Let us compute the bosonized action S[n] l|A.3|) for this time-dependent configuration. The action (|A.3|I is nonlo- 
cal. However, for slowly changing transformations R(t) one can use the "adiabatic approximation" leading to the 
Lagrangian 

L Iot = ^Iab^a^b + W a tt a . (F.2) 
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Here fl a is the angular velocity 

fl = n a t a = -iR^R (F.3) 
and t a are generators of the group G. The parameters of the Lagrangian L Iot are given by equations 

I ab = N c E (i\t a \m)(m\t b \i) + (i\t b \m) {m\t a \i) ^ ^ 

i: occ rn: nonocc 



-N c J2(Ata\i) (F.5) 



w a 

where the summation runs over occupied (?) and nonoccupied (to) eigenstates (|A.6|) of the single-particle Hartree 
Hamiltonian. 

It should be stressed that this "derivation" of the effective Lagrangian l|F.2(l is inspired by the large- N c limit but 
not equivalent to the systematic 1 /iV c -expansion. 

APPENDIX G: MODEL WITH SU(2) FLAVOR GROUP 

The main part of this paper is devoted to the toy model (|45(l with the flavor symmetry SU(3). In this appendix 
we consider the SU (2) analog of this model. Let us take the same Hamiltonian (|37[) 



H toy 



1 2 2 / N c 

/l,/2,/J,/2 = l S 1 ,S 2 ,S' 1 ,S' 2 = 1 \C' = 1 



XV H 2 s',J[ S ' l \hs 2 f 1 S 1 (E^/2 S 2C^/l«cj (G.l) 

but now the indices f' k run only over the values 1, 2. The coefficients Vp s t f s ' |/ 2 s 2 /isi are determined by Eq. i|40[l . 
The diagonalization of this Hamiltonian is quite similar to the case of the SU(3) model. In the baryonic sector we 
have the states with spin J equal to isospin T 

N N N 

J = T = ^, -^-1,-^-2,... (G.2) 

2 2 2 ' 

The spectrum of energies of these states is 

~1 



-N 2 C +2J{J+1) 



(G.3) 



One can describe this spectrum in terms of the 1 /A^-expansion. The Hartree equations corresponding to the I /N c - 
expansion make a subset of equations that we had in the SU(3) case. The single-particle energies £1,62 are given by 
the same equations l|S2*)) . <|%51) . 

The leading 0(N C ) part of the energy coincides with the SU{3) expression l|65|) because it is determined by the 
same expression for the Hartree energy l|79l) . The 0(N C T 1 ) "rotational" J(J+ I) part of the energy also coincides with 
the SU(3) case since it comes from the RPA matrix S. The nontrivial part of the SU(3) matrix S coincides with the 
SU(2) matrix 5*. Therefore the calculation of the moment of inertia, performed in Appendix El for the SU(3) case, is 
also valid for SU(2). The 0(N°) contribution to the energy is given by equation 1|85|) where we must use the levels 
l|HSt : £i = £2, £m = £ i- Since the level S\ has degeneracy 3, we obtain 

E « = -\ E E( £ ™- £ >) =-|(ei-Ea) = §6a2, (G.4) 

m:nonocc i:occ 



which agrees with the 0(N®) part of the exact solution 1JG.3J1 . 
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APPENDIX H: EXTENDED MONOPOLE MODEL 

In this appendix we consider an extended version of the Guadagnini model (|107f) and show that in this extended 
model the picture of "exotic" states relying on the rigid-rotator approximation contradicts the results based on the 
systematic semiclassical analysis corresponding to the large- N c limit. We modify the model I|1U7|) by adding a second 
particle which interacts with the first one via a quadratic pair potential: 

L = ^ | ^ - e(* a • A( Xl )) - y(| Xl |) - Ik (x a - x 2 ) 2 . (H.l) 

Note that this second particle does not interact with the magnetic field and with the potential V. 

Similar to Eqs. (|108|l . (|109|) . we assume the standard semiclassical behavior inspired by the large- N c limit 

Mfc ,e 5 ,y,fc~iV c . (H.2) 

Let us assume that potential V has a minimum at some point R ^ 0: 

V'(R)=0. (H.3) 

Then the static solutions of the equations of motion are described by the conditions 

x (o) =x ( 0)) |x (o )| = jRi (H4) 

Using the rotational invariance we can direct Xj = Xg along the third axis. Now let us write the linearized equations 
of motion for small deviations from the static solution iJx/j (t) = Xk {t) — x[°' . The transverse part (orthogonal to the 
third axis) does not mix with the longitudinal part. It is convenient to use complex coordinates Zk in the transverse 
plane instead of <5x^ 

SxtW^Zke™. (H.5) 
Then the linearized equations of motion take the form 

— /iiw 2 zi + houjzi + k(zi — Z2) = , (H.6) 

— [i2U 2 Z2 — k(zi — z 2 ) = , (H.7) 

where 

h = % . (H.8) 



Equations of motion J 11.6(1 . I(H.7|I lead to the following equation for the frequency 



detl ^ 2 ~^-k k 



k /i2W — k 

Apart from the trivial solution lo = one has three roots of the cubic equation 

u) = k (fiiuj — ho) 1 + (fi2^) 1 ■ (H.10) 

These three roots determine the spectrum of the harmonic excitations. 

Let us show that this spectrum of harmonic excitations differs from the spectrum of the rigid-rotator Hamiltonian. 
Indeed, in the naive rigid-rotator approximation the solution 1|H.4|) leads by analogy with Eq. to the moment 

of inertia of the rigid rotator 

/rot = (Ml + M2) R 2 ■ (H.ll) 

If one inserts this moment of inertia into the Hamiltonian l|112ll and uses the quantization condition l|113() , then similar 
to Eq. H118JI one arrives at the spectrum of harmonic excitations 

E £ = -^-U+^\+0(l/N c ), £ = 0,1,2,... (H.12) 
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corresponding to the frequency 



W rot = — = ; • (H.13) 

Jrot Ml + 



But this frequency is different from the three roots of the cubic equation i|H.10|) . Only in the limit k — > oo (restoring 
the rigidity of the system) one of the roots of equation (|H.10() asymptotically approaches the frequency (|H.13|I . But 
in a general case the rigid-rotator approximation has nothing to do with the true spectrum of harmonic excitations 
which is determined by equations (|102|1 . (|105|l . 
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